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With an aim to argue for the truly collisionless nature of cold dark matter between epochs of 
equality and recombination, we assume a model, wherein strongly coupled baryon-radiation plasma 
ejects out of small regions of concentrated cold dark matter without losing its equilibrium. We use 
the Meschersky equation to describe the dynamics of cold dark matter in the presence of varying 
mass of strongly coupled baryon-radiation plasma. Based on this model, we discuss the growth of 
perturbations in cold dark matter both in the Jeans theory and in the expanding universe using 
Newton's theory. We see the effect of the perturbations in the cold dark matter potential on 
the cosmic microwave background anisotropy that originated at redshifts between equality and 
recombination i.e. 1100 < z < z eq . Also we obtain an expression for the Sachs-Wolfe effect, i.e. 
the CMB temperature anisotropy at decoupling in terms of the perturbations in cold dark matter 
£SJ potential. We obtain similar solutions both in the static and in the expanding universe, for epochs 

of recombination. From this, we infer about the time scale when the dark energy starts to dominate. 
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O I. INTRODUCTION 

u 

Flat cosmological models with a mixture of ordinary baryonic matter, cold dark matter, and cosmological con- 
stant(or quintessence) and a nearly scale-invariant, adiabatic spectrum of density fluctuations are consistent with 
standard inflationary cosmology. They provide an excellent fit to current observations on large scales(>> 1 Mpc). 
Currently, the constitution of the universe is 4% baryons, 23% dark matter and 73% dark energy [TH5]. 

In the standard hot Big Bang model, the universe is initially hot and the energy density is dominated by radiation. 
The transition to matter domination occurs at z ~ 10 4 . In the epochs after equality and before recombination, the 
i— i universe remains hot enough. Thus the gas is ionized, and the electron-photon scattering effectively couples the matter 
and radiation [7J. At z « 1200, the temperature drops below s» 3300 K. The protons and electrons now recombine 
to form neutral hydrogen and neutral helium. This event is usually known as recombination [81-110] . The photons 
then decouple and travel freely. These photons which keep on travelling till present times are observed as the cosmic 
microwave background(CMB). The cold dark matter theory including cosmic inflation is the basis of standard modern 
cosmology. This is favoured by the CMB data and the large scale structure data [TTj, [T^] . The CDM model is based 
on the assumption that the mass of the universe now is dominated by dark matter, which is non-baryonic |13[ 114) . 
Also it acts like a gas of massive, weakly interacting(collisionless)particles [15] ■ They have negligibly small primeval 
I velocty dispersion. Also they are electromagnetically neutral [121 HZ] • 
CNJ The word Cold here means that the ratio ■£ < <j>, the gravitation potential, where T and M represent the temperature 

and mass of the dark matter particle. There is remarkably good agreement between standard CDM models and the 
observed power spectrum of Lyman a observers [18j . This rules out the warm dark matter candidates. The CDM 
model predicts the power spectrum of the angular distribution of the temperature of the 3K cosmic microwave 
background radiation and the flat Friedmann model. A low density CDM model with a density parameter of around 
0. 3 to 0. 4 with cosmological constant actually matches all available data fairly well [HI [20]. The stable CDM 
paradigms predict all structure formation [2TJ [21] . 

The existence of clusters (< 50 mpc) and groups of galaxies suggests that the galaxy formation is due to the 
gravitational instability of a spatially homogenous and isotropic expanding universe. The Perturbations of such a 
model have been first investigated by [231 El] . Then within the framework of general relativity by [25l [26] . Also it was 
studied in a newtonian model [3"2"]. The results obtained from relativistic theory were similar to the standard Jeans 
theory. In the early universe, radiation fixes the expansion rate. This is due to low density and self-gravity of dark 
matter [281 US] • After equality, the main contribution to the gravitational potential is due to the cold dark matter. 
In the metric of the perturbed FRW universe, the main contribution to the gravitational potential comes from an 
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imperfect fluid, i. e. the dark matter. We still consider that the difference tp — <fi is suppressed compared to t/j, atleast 
by the ratio of the photon mean free path to the perturbation scale. Here i\) and <f> are the newtonian potential and 
the perturbations to the spatial curvature in a conformal Newtonian gauge respectively (see eq. 48 sec[n|. We neglect 
the contribution of the baryon density in realistic models, where baryon contributes only a small fraction of the total 
matter density. At ?/ > i] eq , the gravitational potential is mainly due to the perturbations in the cold dark matter. 
The potential of cold dark matter is time-independent both for long wavelength and short wavelength perturbations. 
This is because it is highly non-relativistic. There is a strong coupling between the baryon and the radiation in the 
plasma. So, we treat it as a single perfect fluid for low baryon densities. We neglect the non-diagonal components 
in the energy- momentum tensor of dark matter. This is because for very small photon mean free paths ip is same 
as <f>. Therefore, we treat it as a single perfect fluid for many epochs after equality upto the recombination, ft is 
only after recombination, that the baryon density starts increasing. This is after the primordial nucleosynthesis of 
hydrogen and helium is complete. The radiation after decoupling from matter at very late recombination epochs 
evolves separately. The CMB anisotropy of epochs between equality and recombination gives important information 
about the perurbations in the cold dark matter i. c. the modes that enter the horizon before recombination. 

In our model, the strongly coupled baryon-radiation plasma ejects out of a concentrated region of cold dark matter. 
Based on this model, we study the dynamics of the cold dark matter after equality and upto many epochs near 
recombination. Note that this sort of flow assumes that the baryon -radiation plasma ejects out without losing its 
equilibrium in the presence of cold dark matter. This is possible only if we assume that the dark matter is truly 
collisionless. This assumption can also accomodate a self-interacting dark matter at very small scales, which can 
transfer energy and momentum to the outer core |30j . The strongly-coupled baryon -radiation plasma is only under 
the influence of the potential of cold dark matter at these epochs. We assume that the number of photons per Cold 
dark matter is initially spatially uniform on supercurvature scales (wavelengths greater than H^ 1 ). During this time, 
the matter and radiation densities vary in space. In other words we consider adiabatic perturbations. As the universe 
expands, the inhomogeneity scale becomes smaller than the curvature scale. Thus the components move with respect 
to one another and the entropy of photon per cold dark matter particle varies spatially. In contrast the entropy per 
baryon remains spatially uniform on all scales. This is until the baryons decouple from radiation. 

We use Meshchersky equation to study the dyanmics of the strongly coupled baryon- radiation plasma after the cold 
dark matter starts to dominate. This occurs after epochs of equality. In this model, we imagine a flow of the baryon 
- radiation plasma across regions of highly concentrated non -relativistic cold dark matter. The baryon -radiation 
plasma can be assumed to have an ejection velocity with respect to the cold dark matter after equality. We do not 
disturb the equilibrium of the ejecting baryon-radiation plasma, even in the presence of cold dark matter. This we do 
to argue in favour of truly collisionless nature of cold dark matter. We ensure this by a specific assumption(see eq. 27 
sec. pB. This is the time, when the radiation densities rapidly start to decrease. This is due to its scaling proportional 
to a - , which is a -3 for the cold dark matter density. On the basis of this model, we study the adiabatic perturbations 
between epochs of equality and recombination. This we do, both in the Jean's theory and in an expanding universe 
in Newtonian theory. We then find out the effect of perturbations in cold dark matter potential on the anisotropy 
in temperature of radiation at these epochs. This we do for modes greater than the curvature scales. These are the 
modes which enter horizon before recombination. 

The paper proceeds as follows. In section 2, we discuss the Dynamics of Cold dark matter in Baryon-radiation 
plasma using the Meshchersky equation in the Jeans theory. Then we discuss the adiabatic perturbations in the scope 
of this model. We also deduce an equation which represents how the anisotropy in the temperature of radiation is 
affected by the perturbations in the cold dark matter in the epochs between equality and recombination. We derive 
an expression for the Sachs- Wolfe effect. In Section 3, in the scope of our model, we discuss the dynamics of the 
above mentioned scenario in the expanding universe in the Newtonian theory. We show how the anisotropy in the 
temperature of radiation is affected by the perturbations in Cold dark matter potential, in an expanding universe. 
Also we evaluate an expression for the Sachs- Wolfe effect in the expanding universe. We then write the equation for 
the evolution of cold dark matter perturbations with time. 

II. GRAVITATIONAL INSTABILITY: JEANS THEORY REPRESENTATION 

In Jeans theory we consider a static, non-expanding universe. Also we assume a homogeneous, isotropic background 
with constant time -independent matter density |23l I24j . This assumption is in obvious contradiction to the hydrody- 
namical equations. In fact, the energy density remains unchanged only if the matter is at rest and the gravitational 
force, F cx \j<j) vanishes. This inconsistency can in principle be avoided if we consider a static Einstein universe, 
where the gravitational force of the matter is compensated by the antigravitational force of an appropriately chosen 
cosmological constant. We consider the fundamental equation of dynamics of a mass point with variable mass. This 
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is also referred to as the Meshchersky equation. 



dv „ dm _ 

m dt =F+ ltt U W 



surrounding bodies. The last term ^rrU is referred to as the reactive force. This force appears as a result of the action 
that the added or separated mass exerts on a given body. If mass is added , then the > coincides with the 



It should be pointed out that in an inertial frame, F is interpreted as the force of interaction of a given body with 

dm 
dt 

Tiri-cr An o rri t rnn hi"iri t r It yti aco 10 orlnnrl rn nn tnn 

dt 

vector u. If mass is separated, ^ < 0, and the vector i? is oppositely directed to the vector u. 

We consider a fixed volume element AV in Euler (non-co- moving) co-ordinates x. After equality, when the cold 
dark matter starts to dominate in small regions of space, we write: 

A mS = -AM d ■ Web - V Pbl ■ AV - (2) 
at dt 

Here the first term on L. H. S represents the acceleration in the mass of cold dark matter of mass AMd and Vd is the 
velocity of dark matter element. The first term on the R. H. S. represents the gravitational force on the cold dark 
matter. The second term on right is for the force due to the pressure of the baryon -radiation plasma .31j. The last 
term on R. H. S. is the reactive force on the cold dark matter due to the ejection of baryon - radiation plasma from 
regions dominated by cold dark matter. This last force arises only due to the model that we assume here. Here 
is the pressure of the baryon-radiation plasma and u is the ejection velocity of the Baryon - radiation plasma with 
respect to the concentrated region of cold dark matter. Note that in (eq. 2), we do not take the pressure of cold 
dark matter into account. This is because, the cold dark matter is highly non - relativistic. Therefore, we neglect its 
pressure. We assume that after equality, the strongly coupled Baryon - radiation plasma starts to rapidly decouple 
from the matter. The matter at these epochs is predominantly Cold Dark matter. 

After Equality, in regions dominated by highly non - relativistic Cold dark matter, we write the continuity equation 
for the ejection of baryon- radiation plasma complex : 



dMb-y f deb-f 



dt ./av dt 



dV (3) 



where Mf, 7 is the ejecting mass of baryon-radiation plasma and e& 7 is the energy density of baryon radiation plasma in 
the concentrated region of cold dark matter, from where it is ejecting out. In this model, we assume that the ejection 
of baryon -radiation plasma out of a concentrated region of heavy Cold dark matter does not disturb the equilibirium 
of the dark matter. This we can assume only because of the truly collisionless nature of Cold dark matter. This is 
the fundamental assumtion of this model. We neglect the flux of the Cold dark matter out of a region of volume AV 
in the time that the baryon -radiation plasma flows out of this region. The rate of flow is entirely determined by the 
flux of the baryon - radiation plasma and we write : 



We write (Eq. 2) as: 



e d AV^ = -e d AVV4> - V Pbl AV - ^ AVu (5) 



where e d is the energy density of cold dark matter, and <j) is the gravitational potential of cold dark matter. 

r dvd _ _ _ , „ dsb-y _ ,„, 
£ d[-Qj- + v d ■ Vu d J = -£d^4> - Vp b7 gj-U (6) 

Now using the Jeans theory, we introduce small perturbations about the equilibrium values of variables, W. Bonnor 32J: 

£ d (x, t) = E do + 5e d {$, t) (7) 

Vd{x, t) = v2o + Sv d (x, t) = 5vd(x, t) (8) 



<p(x, t) — O + S(j>(x, t) 



(9) 
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u(x,t) = u D + 5u(x,t) (10) 
where v do « c, the speed of light, u a ^ and Se d « e do 

Ptry(x,t) = p bl {e blo + 5e bl , S + SS) = p blo + 5p bl (x,t) (11) 

S(x,t) = S + 6S(x,t) (12) 
where S is the entropy of cold dark matter element. Also we write: 

dpb-y = c 2 5e b7 + aSS (13) 
where c 2 is the speed of sound. Neglecting Dissipation, we write: 

dS dS _ - 

dt = m +Vd ' vs (14) 

V 2 = AnGs d (15) 

From (eq. 4) we write: 

°^+V-e hl u = Q (16) 

We substitute the values from (cq. 7)to(cq. 15)in(eq. 16)and(eq. 6) to get: 

dSe b -y 



dt 



+ e blo V -Su + V ■ (SshyUo) = (17) 



£ do^^~ = -5e d V(t> - V{c 2 s 5e bl + a8S bl ) - ^j^-Su- ^^-u (18) 

V 2 S(j) = A7rG5s d (19) 
dSS dSS ,._ _. _ „ 

— = — + m .V)S o = (20) 



We now take the divergence of (eq. 18)to get: 

£do — ~Q t — = -(V • 5e d )Vcj) - Se d {ATrGe do ) - c 2 s \7 2 Se bl 

2 - d(e bjo Su) dSu - d{5e bl u ) du a 

_ ctV 5S-V- e blo — } V • {^^- - 5e bl — } (21) 

Now using (eq. 17) and (eq. 19) we get : 

- c 2 s V 2 6s b , 4nGs do Se d (V • fe d )V0 o - e j^ ^ V " d) + V • {e b J-^ + Ss^} = aV 2 5S (22) 

In a static universe the total energy remains constant. So for small inhomogeneities we write: 

<5e b7 = -6s d (23) 

Thus (eq. 22) becomes: 

- c 2 8 V 2 8e bl + 4nGe do Se bl + (V • & h7 )V0 o - e do — + V • + fc <>7^} = oV 2 S{x) (24) 
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For strongly -coupled baryon - radiation plasma before recombination, with low baryon densities we write : 



Thus we get: 



£b 7 + Pbf — £b + g£ 7 — g^2 £ f (^) 



4 

6s bj + 5p bl = (26) 



The baryon - radiation plasma is only affected by the gravitational potential of the cold dark matter in these epochs. 
We argue that the baryon-radiation plasma ejecting out of concentrated regions of cold dark matter is always in 
equilibrium, even in the presence of cold dark matter. This is possible only if the cold dark matter is truly collisionless. 
So for epochs between equality and recombination, we write: 

dPb-y _ £fc 7 + Pb-y / 07 n 

dT ~ T {V 

The above equation is of profound importance in this model. This is because, it ensures that the equilibrium of the 
ejecting baryon -radiation plasma is not disturbed, even when it flows out of concentrated regions of cold dark matter. 
This we can write, because the distribution function of pressure of the baryon-radiation plasma depends only on ^. 
We assume here that the chemical potential is much smaller than the temperature. So from small perturbations about 
the equilibrium values of variables p bl andet 7 in the (eq. 27), about a fixed temperature T, we get, 



Spb-y 4(5e 7 
~5T = 9^f 



(28) 



Thus we write: 



4&^ ra ST. 

(29) 

Here ^ = O, which represents the anisotropy in temperature of radiation for small baryon densities after equality 
For small perturbations in cold dark matter, using (eq. 29) in (eq. 24), we write: 



^ 2 fe 7 {l 9} _ 4 V 2 {(M1 _ e}} + v 2 o Afe 7{1 _ e} + Ay • {5e.il 6)}V0 O 

-eJ^- + V • {ebjf} + • {fey (l - = *V>5S( X ) (30) 

For small perturbations in Cold dark matter, we neglect the term V • 5v*d- This is due to cold dark matter being 
highly non -relativistic. Also we neglect V • ^f- . This is due to the fact, that in this model, we assume that the ejecting 
velocity of baryon - radiation plasma does not suffer interactions and collisions. This is because of dark matter being 
truly collisionless. So, we further write: 

4 3 2 fe 7 {l-9} 4 V 2 {&7{1 _ e}}+47rGe<to 4 _ Q} 



9 C 2 dt 2 9 L 7L JJ 9c 2 

4 - 4 -. 81T 

+ g^ V ' - ©)Wo + • {fe 7 (l - 6)^} = aV 2 SS(x) (31) 

In the above equation, we see the contribution of divergence in the anisotropy of radiation at epochs between equality 
and recombination. At epochs near recombination, the strongly coupled baryon-radiation plasma starts rapidly to 
decouple from the cold dark matter. During these epochs, therefore the actual velocity of the baryon -radiation 
plasma does not change appreciably due to the gravitational force of cold dark matter. During these epochs, the rate 
of increase of ejection velocity is only due to gravitation potential of cold dark matter. There is no force due to the 
pressure of the baryon -radiation plasma. Therefore we write: 
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We can argue that due to high baryon densities, the pressure of the baryon-radiation plasma vanishes at these 
epochs. But the baryon densities can never be very high. This is because it will hinder hydrogen nucleosynthesis 
after recombination. Thus it has to be assumed that the pressure of the baryon - radiation plasma vanishes not due 
to very high baryon densities, but due to a strange form of energy(dark cnergy)with negative pressure, which has 
started to dominate near recombination. Therefore, the acceleration in the ejection velocity will be primarily due to 
the gravitational forces So, we write: 

u = v^f - v2o (33) 
u = Vf^ - v2o (34) 

vjy = (35) 

The above (eq. 35) shows that the actual velocity of the baryon-radiation plasma stops to increase at epochs, where 
it is almost about to decouple from cold dark matter. Therefore, at these epochs the accleration of cold dark matter 
reverses its sign and we write: 

v2o = V<^ (36) 

With these Assumptions the (eq. 31) reduces to: 

^{fc 7 l - G}dt 2 - c 2 s V 2 {Se 7 (l - &)} + 4nGe do {5s, ( (l - 6)} = -±*V 2 8S(x) (37) 
Considering adiabatic perturbations: 



We use: 



Then we write: 



The (eq. 37) now reduces to: 



SS = (38) 



d 3 k 

*:,0 - (-))(-■/)= / 5e 7k (t)(l-Q) k (t)e ik -*j^ (39) 



6s 7k (t)(l-e) k (t)=y k (t) (40) 

Vk(t) + k 2 c 2 s y k {t) + ATrGe doVk {t) = (41) 
or 



y k (t) + (k 2 c 2 s + 4irGe do )y k (t) = (42) 
The above equation has two solutions:^ oc e ±IW 0, where 



Here the Jeans length is : 



u(k) - ^JWcl + AnGe do = Jk*c* - ^Ge blo (l - — ) (43) 



i3 (44) 



Thus uj(k) is real for A < Xj or k > kj, where 



kj _ (45) 



C s 
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when 

k 2 c\ > ^Ge blo {\ - — } (46) 

Where e is the total equilibrium value of energy density and £b 70 is the equilibrium energy density of baryon-radiation 
plasma, at these epochs. We interpret from above, that all the modes near recombination are real. This is because the 
value of starts increasing rapidly near recombination. This is due to the decoupling of baryon - radiation plasma 

from cold dark matter. This explains the maximum frequency of fluctuations in the CMB temperature anisotropy at 
red -shifts of recombination. For A < Xj the solutions are: 

yk(x, t) oc sin(u;£ + kx + a) (47) 

For A < r 7 < Xj, where r 7 stands for the mean free path of free- streaming photons after equality, free-streaming 
becomes important. Free-streaming refers to the propagation of photons without scattering. We write the equation 
for free-streaming photons in the conformal Newtonian gauge [33] . The gravitational potential is primarily due to the 
cold dark matter. We argue that the photons can still be described by the equilibrium distribution functions. This is 
because the cold Dark matter is collisionless. Therefore, their mutual interactions will not disturb the equilibrium of 
the ejecting baryon-radiation plasma. Also, we treat the baryon-radiation plasma as free-streaming photons for low 
baryon-densities. We use the metric below: 

ds 2 = -1 - 2?P(x, t)dt 2 + a 2 5 l3 (l + 2<j)(x, t))dx l dx j (48) 

where ip corresponds to the Newtonian potential and <p is the perturbation to the spatial curvature. The equation for 
free- streaming photons is: 

I^ = _#_^_^^ f49 ) 
pdt dt a dx i v ; 

ip = <p (50) 

We can make the above assumption, for epochs between equality and recombination. This is because after Equality, 
and before recombination, the dominant contribution to the potential is due to cold dark matter. So, we neglect the 
non-diagonal components of the energy-momentum tensor of the baryon-radiation plasma. Recall that the difference 
ip — <p is suppressed compared to ip, atleast by the ratio of the photon mean free path to the perturbation scale. 
Thus we assume it to be a single perfect fluid. In a static Einstein universe, we neglect the terms due to the Hubble 
parameter H . At scales of the order of mean- free path of the free-streaming phtons in the baryon-radiation plasma, 
i. e. of the order of r 7 , we neglect the spatial inhomogeneities We therefore write: 

1 dp _ d(p 

pdt~ dt (5ij 



Using the assumption of (eq. 50), we write: 



1 d P d ^ /CO) 

p dt - dt [bA) 



Recall that the Newtonian potential ip corresponds to the dark matter potential <fid - This <t>d remains constant 
between epochs of equality and recombination. Thus we write: 

\ogp = -ip (53) 

\ogp = -4> do ^P = e~ 0do Sp = -e'^Sfa (54) 

The momentum per unit volume of radiation at epochs of decoupling is equivalent to radiation pressure. Also, if the 
photons, just at the epochs of decoupling are in thermal equilibrium, and that it has same average energy associated 
with each independent degree of freedom, we write: 

§P 1 = ^Sp (55) 
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where 5p 1 is the increment in radiation pressure associated with each independent degree of freedom at decoupling. 
Since for a unit volume decoupling photons' energy density is same as radiation pressure per degree of freedom, we 
write: 

Se y = ~e-*"6<l> d (56) 

We use (eq. 40 and (eq. 47) to write: 

& 7fc (l - Q)(x,t) = Asin(ojt + kx + a) (57) 

where 



w= Jfc 2 c2-47rG£ 670 (l- — ) (58) 



Therefore, for epochs of decoupling, we write: 

-e-**°66 d 



{1 - 6} = Asm(ut + kx + a) (59) 



3 

or 



From the above equation we can write: 



Therefore, we write: 



3 e <pdo 

{0 - 1} = — A sin(wi + kx + a) (60) 

0(f>d 



^ = {e + 4> d {i + ef} (6i) 



i<Pd 

3 



(62) 



The above result is the same as that predicted by Sachs- Wolfe. Three types of effects(due to fluctuations in density, 
velocities and potential)simultaneously contribute to the CMB temperature anisotropy. The fluctuations that matter 
at scales beyond 1° are those in the gravitational potential 8(f) d (Sachs- Wolfe effect) [34], 

The above equation represents the anisotropy in the CMB temperature, for epochs near recombination, for regions 
where sufficient primordial helium synthesis takes place, even before recombination, or the dark energy has started to 
dominate. Recall that in assumption of(eq. 32), we argue that the force due to pressure of baryon-radiation plasma 
vanishes near recombination epochs, and that it can be only due to the fact that the dark energy with negative 
pressure has started to dominate. The value of (f>d is constant. This is because the potential of cold dark matter 
remains constant for many epochs between equality and recombination. Thus in a static universe, with this model, 
we see that the dominant component in CMB temperature anisotropy fluctuations is near recombination. This is 
because, at epochs near recombination, the perturbations in the cold dark matter potential are very small. This is 
because in this model, the decoupling of the baryon - radiation plasma from concentrated regions of cold dark matter, 
near epochs of recombination is almost near completion. 

We now write the CMB fluctuations for supercurvature modes i. e modes with A >> A,/, and for regions where 
sufficient primordial helium synthesis takes place, even before recombination, i. e. the baryon-densities are high. We 
neglect the effect of gravity at these epochs. We can do so because at late recombination epochs, the pressure of 
baryon-radiation plasma is low, and the negative pressure of a strange form of energy (dark energy), which starts to 
dominate at these epochs, cancels the effect of forces due to low pressure baryon-radiation plasma and that of gravity, 
we therefore write: 

Vk oc e ±lkc ° l (63) 

Also for A >> Aj, free-streaming of photons is no longer relevant. Therefore the scattering of photons will dilute the 
anisotropy to a large extent. For late recombination epochs, when the radiation has almost decoupled from matter, 
we write: 

c* = - =*> y k oc e ± 7f (64) 
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We have neglected gravity here again, because at these epochs, the effect of dark energy, which had started to 
dominate from some earlier epochs is to cancel the forces due to gravity and the pressure of baryon-radiation plasma. 
This is possible even in regions of lower baryon densities, which has higher pressure than the regions of higher baryon 
densities. This is only because the dark energy with negative pressure had been dominating from some earlier re- 
combination epochs. The above equation shows that the frequency of fluctuations in CMB temperature anisotropy 
spectrum of supercurvature modes (modes which enter the horizon very early near recombination epochs), at late 
recombination, remain constant till today. This is valid both for regions of low or higher baryon densities (where suffi- 
cient primordial helium nucleosynthesis takes place before recombination). This is because at very late recombination 
epochs when the radiation has almost decoupled fully from matter, the speed of sound approaches a constant value 
of |. The effect of scattering of photons will dilute this anisotropy in supercurvature modes. So, it is of not much 
cosmological significance. 



III. INSTABILITY IN EXPANDING UNIVERSE: NEWTONIAN THEORY 

Using our model (see Sec. [I]), we study the same scenario, i. e. of cold dark matter in the presence of strongly 
coupled baryon -radiation plasma at epochs between equality and recombination. Here we use Newtonian theory of 
expanding universe |32j . We treat the dynamics between cold dark matter and the baryon - radiation plasma, again 
in the framework of the Mcshchersky equation. We assume that the strongly coupled baryon -radiation plasma is in 
the presence of a gravitational potential. This potential, is only due to cold dark matter, at epochs after equality 
and before recombination. We neglect the non-diagonal components of the energy-momentum tensor of the cold dark 
matter. This is because, the difference ip — <p is suppressed as compared to tp, atleast by the ratio of the mean free 
path to the perturbation scale. We consider dark matter as highly non -relativistic fluid compared to the baryon - 
radiation plasma. We consider epochs when the Dark matter has already started clustering. We neglect the rate of 
flow of cold dark matter out of a given region of space, as compared to the baryon-radiation plasma. Therefore we 
write: 

^+W-(e bl u)=0 (65) 

where u is the relative velocity of the baryon - radiation plasma with respect to the cold dark matter. We write 
the Mcshchersky equation for cold dark matter with baryon radiation plasma, ejecting out of concentrated regions of 
space dominated by cold dark matter: 

£do{^^ + {vdo ■ V)^} = -£dcN(t>o ~ Vp h7 - -^p u "o (66) 
In an expanding flat, Isotropic and homogenous universe, we write: 

«b* 7 = v blo (t)]e bl = e b70 {t);u = H(t).x;e do = £do{t)\v2o = v do (t) (67) 

Therefore we write: 



de blo {t) 



+ e blo (t)V ■ u = (68) 



We first discuss for scales >> H^ 1 i. e the curvature scale. In this case, we neglect the velocities of Cold dark 
matter particles. This is because, there is insufficient time to move highly non -relativistic cold dark matter upto 
distances greater than the Hubble scale. Therefore the entropy per cold dark matter is conserved on supercurvature 
scales ((krj > 1)). So we write: 



de blo (t) 
dt 

The (eq. 63) gives: 



+ 3e blo (t)H(t) = (69) 



Bf 

SdoVcbo + Vp fc7 + -t^u = (70) 

We take the divergence of the above equation to write: 

6H 2 e byo = AnGe 2 do (71) 
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To get the above equation, for near recombination epochs, we use: 

V Vy = (72) 

We use the above assumption because for epochs near recombination, the baryon-radiation plasma starts to decou- 
ple fast. Therefore, only the gravitational force determines the acceleration of dark matter. Now for Adiabatic 
perturbations, we write: 

SS = (73) 
Also, we write the following perturbations for other variables: 

£d = £d = £do + Se d (x, t);u = u + Su(x, t);<p = <p + 5(f>;p bl = p blo + 8p bl = p Q + c 2 s 8e bl (74) 
Where the variables have their usual meanings (see Sec. [n]). We use (eq. 71) in (eq. 66) and (eq. 68) to write : 

^ = £b 70 V • Sit + V ■ {8e bl u } = (75) 
and 

e do V5cj) + 5e d V4> + V<% 70 + %^tZ + ^g^o = (76) 

at at 



We use the Langragian co-ordinates and write: 



{§~ t h = {^h ~ to ■ V, (77) 



where 

u a = H(t)q; x = aq (78) 

We write: 

1 Se 

V x = -V q ;S=- (79) 

d Eq 

where S is the fractional amplitude of perturbations. We write (eq. 75) in the Co-moving coordinates, using (eq. 65) 
and (eq. 77) to get: 

<5 b7 + = (80) 

a 

also we write the (eq. 76) in the co-moving coordinates to get: 

^ do ^do ^ ' -* 

— V(S0 H -S d V(t) + —e^oVShy + e bjo Su + e blo S b7 + e blo 5 bl -e bjo \7 ■ S by — (81) 

a a a a 



We write (eq. 69)in co -moving coordinates: 



from (eq. 71), we get: 



ie blo {t)H{t) 

£byo — (O^) 



We take the divergence of (eq. 76), and use (eqn's. 68-77-79) and (eq. 80) to get: 

+ 3aH5 bl + { t^}V% 7 + V^ 7 - V.<5 b7 { g(2 " + 3) } = (84) 
a a a 
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We Neglect the fourth term i. e. V.#6 7 . This is because, it is proportional to v,Vl5 " , which is very small. This is 
due to very small divergence in the perturbations to the ejection velocity of baryon -radiation plasma. With this 
assumption, again we argue that the cold dark matter is truly collisionless. So, in this model, it does not disturb the 
equilibrium of the ejecting baryon-radiation plasma. From (eq. 71), we get: 

8irGe 2 do 6 d = 6H 2 e 2 „ (0 5 bl (85) 

Using the above, we write (eq. 84) as: 



Therefore, neglecting the term containing V.<5f, 7 , we write: 



We use: 



Thus, we write the equation for the evolution of fractional amplitudes of perturbations in the baryon -radiation plasma 
as: 

2 _ 2 r \H 3 F 2 

5 hl = exp{fc 2 (^^) - ^^}t (89) 

For very late recombination epochs, the radiation starts to decouple rapidly from matter. Therefore, baryon densities 
are very low. So we write: 

c 2 s = \ (90) 

where c 2 is the speed of sound. 

The (eq. 89) shows, that for low baryon densities, fractional amplitudes of perturbations in the radiation, which 
originate at late recombination epochs grow at the fastest rate. This is because the value of c 2 is maximum at late 
recombination epochs. It is because of the lowest baryon densities, in the coupled baryon -radiation plasma at these 
epochs. Also we can see that the second term in the exponent in (eq. 89)vanishes for very late recombination epochs, 

as ratio -Sp- — > 0. We thus write the equation for evolution of the fractional amplitude of perturbations in the Cold 
dark matter, which originate after equality, when the clustering of dark matter had already started. 

•? ff2 p 2 2 ,,2 1FT 3 p 2 

= " (91) 

4irG e do 6a A H 8irGa z e do 

The above equation represents the growth of the fractional amplitudes of perturbations in the cold dark matter, 
which originate after equality, for scales >> H . If baryons contribute a significant factor of the total matter 
density, CDM growth rate will be slowed down between equality and the recombination epochs 35 . Also we see 
that the perturbation growth rate will slow with scale factor [35]. The CDM density fluctuations will dominate the 
density perturbations of baryon-radiation plasma |35j . This is because for scales >> H^ 1 , the density perturbations 
of baryon-radiation plasma are washed out by the scattering of photons at scales > r 7 , which is the mean free path of 
photons. The perturbationa in the cold dark matter will cease to grow when the dark energy starts to dominate [37] . 
We can interpret this from (eq. 91). This is because with the growth of dark energy, the value of c 2 will decrease. 
This will then lead to ceasing of growth of CDM perturbations. There is existence of non-linear structures today. This 
implies that the growth of fluctuations must have been driven by non-baryonic dark matter, which was not relativistic 
at recombination. Also, we see that the perturbations at supercurvature scales grow slowly. Recall that these are the 
modes which enter the horizon very early, well before the recombination epochs. The slow growth of such modes is 
because, it is only at late recombination epochs that the second term in the exponent in (eq. 91) will vanish due to 
— > 0. Also, it is only at late recombination epochs, that the value of c 2 reaches its maximum value of |, just 
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before decoupling. The amplitude of the fractional density perturbations in the cold dark matter, in (eq. 91) will be 
maximum when the ratio — > 1. This will occur when 

°*> = 2ttG (92) 

In writing the above equation, we use the result of (cq. 71). So for epochs when the density of baryon-radiation 
plasma is equal to the density of the cold dark matter, (eq. 91) is: 

fc-T-^w'-S* (93) 

Now we discuss the gravitational instability for scales < H -1 . This originates after equality with dark matter 
dominance in the presence of strongly coupled baryon - radiation plasma in an expanding universe. We discuss it in 
the Newtonian theory. The decoupling of the strongly coupled baryon - radiation plasma from the non -rclativistic 
cold dark matter starts after equality. Let us assume that the separation of this plasma from the cold dark matter 
gives a relative velocity of u Q to the baryon - radiation plasma. At small scales with low baryon densities, we neglect 
the contribution of non-diagonal components in the energy- momentum tensor of dark matter. Therefore, we treat 
it as a perfect fluid for many epochs between equality and recombination. We treat the strongly coupled baryon - 
radiation plasma as a perfect fluid for epochs between equality and recombination. This is because of low baryon 
-densities at these epochs. Recall that the baryon densities only starts increasing substantially after recombination. 
This is when the primordial nucleosynthesis of hydrogen and helium will start. However, as an exception in certain 
regions, the primordial nucleosynthesis may start at epochs before recombination. We now conclude that there is 
sufficient time for the Cold dark matter to flow through distances at scales < H~ x . This is because from (eq. 91), we 
see that the perturbations in the cold dark matter, at small scales, grows at the fastest rate. Therefore we write the 
Meshchersky equation as below: 



£do{^7^- + {via ■ y)vdo} = -e do V(t> - Vp bl - ^^-Uo (<)4 ) 

where 

u = Vb 10 - v2o (95) 

At small scales, we assume that in the time that the baryon -radiation plasma flows out of a given region of space, 
the inhomogeneity in Cold dark matter in that time duration is negligible. So we assume that the baryon-radiation 
plasma can flow out of a concentrated region of dark matter without generating the collision terms. This is because 
the cold dark matter is highly non-relativistic and collisionless. Also, because the dark matter has already started to 
cluster. Therefore we write: 

v~b-t = v b70 (t) = H(t)x (96) 

where x is the eulerian cordinate. We write the continuity equation for flow of the strongly coupled baryon - radiation 
plasma as : 

F)f 

+ V.(£6 7 «b 7 ) = (97) 

Here £f, 7 is equal to £b 70 (£) 

£bjo — ~3Hebyo (98) 

We take the divergence of (eq. 94) and neglect the spatial dependence of u . This is due to the collisionless nature of 
dark matter. We also assume that e do is constant for small scales. We thus write the Friedmann equation: 

H + H* = - A * G£do (99) 
3 

For adiabatic perturbations we write from : 

s d (x,t) = £ do + Se d (x, t); v d (x,t) = v2 Q + 5vd(x,t) = 5v d {x,t); 4>(x,t) = 4> a + 84>{x,t); 

u(x, t) = u + Su(x, t);p bl = p b7 o + Cs&t 7 ; SS = (100) 
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where the variables have their usual meanings (see Sec. [TT]) . Then we use the perturbed values of variables in (eq. 
94), to write: 

Sdodvfry - e do 8u + 5e d vC 10 - Se d u + e do (v^y • V)£i£y - e do u ■ V8vb 1 - e do v^ ■ 

+e do u ■ V5u = -e do S754> - 6e d V4> - c 2 s Ve bj - e b \ Su ~ Se b ~/U (101) 

We then use the perturbed values of variables in (eq. 97) to get: 

Se'b-y + £b 7 o(i)V .Svb 7 + Ufc^oV.fe^ = (102) 

We now write (eq. 98) in co-moving coordinates. 

£ ^ o = _^^ (103) 

We also write (eq. 102) in co-moving coordinates to get: 

&' fc7 = -^V. ( 5w b ' 7 (104) 

Using (eq. 95) in (eq. 94) and also using the (eq. 103), (eq. 104), we write (eq. 94) in co-moving coordinates. We 
then take the divergence of the obtained equation and use the results of (eq. 112) and (eq. 117), to write: 

•■ 3 u a ■ 3H 3H 2 eb 7 o , 9H 2 £b- fo 3Hu Q 47rGe& 70 

d b7 + H { H — }<5 b7 - { 1 — H 1 5 h 4vrGe fc70 — }d bl 

a atl a e do ae do a z a 

— — —{av blo - 1 - u \V.5 bl ■ H -{ y }V 5 bl = (105) 



a £do ae do e do a' 



Also we can write: 



_ 



where 9 = represents the anisotropy in the temperature of radiation. The total energy distribution composed of 
the sum of the dark matter and the baryon-radiation plasma is a function of time in an expanding universe. If this 
remains smooth inspite of the inhomogeneities at scales smaller than H —1 , we can write: 

s bl a 4 {t) + e d a 3 {t) = E total (t) (107) 

or for any time 't' , we can write: 

— aSeb-y — Se d (108) 

At epochs near recombination, when the baryon-radiation plasma is decoupled from the cold dark matter, we assume 
that their bulk velocity does not appreciably change in these epochs. Therefore, it remains constant for epochs near 
recombination. The ejection velocity of baryon-radiation plasma is affected only due to the gravitational force of the 
dark matter potential. There is no force due to the pressure of the baryon-radiation plasma. This can be either due 
to low pressure of the plasma, due to higher baryonic densities, or due to some new force of strange energy(dark 
energy), having negative pressure. The baryon densities can never be very high before recombination. This is because 
it would hinder future hydrogen nucleosynthesis after recombination. So, it would be rightly concluded that a new 
form of dark energy having negative pressure starts to dominate at these epochs. So, we write: 

Vb~to = constant => u&* 7 = (109) 



Vdo = V0 O (110) 

The above equation comes from the same argument, which we give in (eq. 36). The actual velocity of the baryon- 
radiation plasma stops to increase, when it is almost about to decouple from dark matter. At these epochs, the 
acceleration of dark matter will reverse its sign. Using the above equation and (eq. 95), we write: 



u = -V0, 



(HI) 
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and in co-moving coordinates, 



£ = -^ (112) 



V.Sii = -V6<j> (113) 

and in co-moving coordinates, 

V.*« = -^ (114) 

or 

V.W = -^ = -^ (115) 
a a 

Using (eq. 103), we can write the above equation as: 

V.Sii = AnGeb-yodb-y (116) 

Also we assume, 

V.Su = (117) 

For scales much smaller than the Jeans lengthA < t 7 < Aj i. e. scales much smaller than the curvature scale i. e 
for which H~ x is dominant, we neglect the terms oc H and write Eq. 105 as: 



Vk -Vk + {4^? £b70 (l - 1 ) - fc 2 e fc7Q ( C ' V ?° U ° ) - ^}y k = (118) 

a a a 2 a 

We use (eq. 39)and ((eq. 40) in writing the above equation. If we now choose: 



2A 



Uo 

a 



{47rGe b70 (l - ±) - k 2 s b , { cl ^°"° ) - ^} = u, 2 (119) 

Then the auxiliary equation is: 

D 2 + 2XD + uj 2 = (120) 

with 

D = -X±V\ 2 -lo 2 (121) 

The only solution which is of physical significance is when A < w, i. e. when the roots of the auxiliary equation are 
imaginary, i. e. : 

D = -A + ia 

2 2x2 

a = lu — A 

y k = A\ {I + (-) 2 }e~ At cos{a - arctan-} (122) 
V a a 

Thus we write: 



<):-,{ \ -(-)} = u/{l + (^) 2 } e - At cos{a- arctan^} (123) 



for A < t 7 < Aj, we can write for free -streaming photons: 
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where tp corresponds to the Newtonian potential in Perturbed FRW universe. Here we assume, 

</> = i\) (125) 

This is because the potential is dominated by the cold dark matter between epochs of equality and recombination. 
We neglect the spatial inhomogcncities in the Newtonian potential at very small scales and terms oc H . Using the 
same argument that we give for (eq. 56. ), we get: 

fc 7 = -Sp = (126) 

Using (cq. 123)and (eq. 126), we write: 



= A\ {I + (-) 2 }e~ At cos{a - arctan -} (127) 
1 a a 

or 



{6-1} = ^r— \/{l + (-) 2 }e" A * cos{a- arctan-} (128) 
oq>d V a a 

This equation is of same form as the one we get for Einstein's static universe in Jeans theory. In the static universe, 
we arrived at this equation by an assumption of (eq. 32). This we reasoned, may be partly possible where baryon 
densities are higher due to early primordial synthesis of Helium before recombination. But the primary reason would 
be that dark energy, with negative pressure starts to dominate from these epochs. The similarity in form of the two 
equations, shows that the epochs, for which the CMB temperature anisotropy derived from the solutions in the static 
Einstein universe and those from the expanding Newtonian match, will be the epochs when the dark energy starts 
to dominate. We see that the amplitude in (eq. 123)is variable, but at epochs when the temperature anisotropy of 
radiation is dominant, which is the late recombination epochs, the term A can be negleted. So we write (eq. 123) as: 

{6-1} = ^— cos a (129) 



Using the result of (eq. 61), we write: 



f (130) 



The equation for fractional amplitudes of perturbations in the cold dark matter that originate at scales smaller 
than the mean free path of the photons is: 



5 d = + (-) 2 }e" A * cos{a - arctan -} (131) 

9e do V a a 



and for late recombination epochs, neglecting A, we can write: 



AAc 2 a 

8 d = -^p^cosa (132) 
y£do 



where 



The above equation shows that the fractional amplitudes of perturbations in cold dark matter at small scales after 
equality keep growing with scale factor, but with a decreasing frequency of oscillations. 
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IV. CONCLUSION 



We conclude that this model, wherein we assume that the baryon -radiation plasma has an ejection velocity, flowing 
out of regions of highly non- relativistic cold dark matter, explains the effect of cold dark matter perturbations on CMB 
temperature anisotropy fluctuations fairly well. The results of (eq. 62) and (eq. 130) represent the Sachs- Wolfe effect. 
It also correctly predicts the growth of fractional amplitudes of perturbations in cold dark matter, which originate 
between equality and recombination epochs. The (eq. 59) and (eq. 127) describe the affect of the perturbations 
in the cold dark matter on the CMB temperature anisotropy fluctuations. They also highlight that the dominant 
component of CMB temperature anistropy fluctuations is at red-shifts of very late recombination epochs z sa 1100. 
This is because the perturabtions in the cold dark matter are minimal at these epochs. The unperturbed potential (f>do 
remains frozen at epochs between equality and recombination. Therefore the fluctuations in the CMB temperature 
anisotropy between epochs of equality and recombination are determined by the perturbations in the potential of 
cold dark matter potential only. The assumptions of (eq. 32) and (eq. Ill), assume the growth of a strange form of 
energy (dark energy), with negative pressure, which cancels the force due to the pressure of baryon-radiation plasma. 

These assumptions also represent the fact that the rate of increase of ejection velocity of the Baryon -radiation 
plasma is only due to the gravitational field of the Cold dark matter. It can be inferred that the force due to its pressure 
can vanish only to some extent due to low pressure in high baryonic density regions. This is due to early primordial 
helium nucleosynthesis in such regions. Since the baryonic densities can never be very high before recombination. 
This is because it will hinder the hydrogen nucleosynthesis which starts after recombination. Therefore, it has to 
be assumed that the primary reason for the vanishing of pressure of baryon-radiation plasma, at near recombination 
epochs, is that the dark energy term with negative pressure starts to dominate at these epochs. The epochs when 
the solutions of the (eq. 59) and (eq. 127) match will give the correct time scales when the Dark energy starts 
to dominate. This work can be taken up in a following paper. The use of (eq. 27) for strongly coupled baryon 
-radiation plasma indicates its equilibrium at all epochs between equality and recombination, even in the presence of 
dark matter. This assumtion can only be made if the dark matter is truly collisionless. The correct explanation of 
pridominant contribution to CMB temperature anisotropy from very late recombination epochs, and the Sachs- Wolfe 
effect, based on the assumptions of (eq. 27), emphasizes the truly collisionless nature of cold dark matter. 

Thus it can be concluded that the use of this model for describing the dynamics of Cold dark matter, in presence 
of strongly coupled Baryon -radiation plasma for epochs between equality and recombination is sufficiently useful. 
More accurate results can be obtained by studying the dynamics of baryon -radiation plasma using energy momentum 
tensor of an imperfect fluid in such a scenario. This will be attempted in a future work. 



Acknowledgments 

We wish to thank Md. Mahfoozul Haque for valuable advice in formatting of the manuscript in Latex. Also, we 
thank Zaheer and Divyendu for valuable discussions and reading the manuscript. 



[1] B. Reid et al., Cosmological constraints from the clustering of the Sloan Digital Sky survey DR7 luminous red galaxies, 
MNRAS, textbf404, 60 (2010). 

[2] W. Percival et al., Baryon accoustic oscillations in the Sloan digital Sky survey data release 7 galaxy sample, MNRAS, 
401, 2148 (2010). 

[3] E. M. Huff et al., Seeing in the Dark ILCosmic shear in the Sloan Digital Sky survey, arXiv:1112. 3143 
[4] J. L. Tinker et al., Cosmological constraints from galaxy clustering and Mass-to-number ratio of galaxy clusters, APJ. 
745, 16 (2012). 

[5] D. Larson et al., Seven year WMAP observations: Power spectrum and W-map derived parameters, APJS. 192, 16 (2011). 
[6] J. Dunkey et al., The atacama cosmology telescope:Cosmological parameters from the 2008 power spectrum, Apj. 739, 52 
(2011). 

[7] P. J. E. Peebles, J. T. Yu, The Astrophysical Journal. 162:815-836, (1970). 

[8] P. J. E. Peebles, The Astrophysical Journal, 153, (1968). 

[9] Y. zel'dovich, K. V. Sunyaev, Sov. Phys. -JETP, 28:146 (1969). 
[10] S. Seager, D. D, Sasselov, D. Scott, Ap. J. sup. 128:407430 (2000). 
[11] P. J. E. Peebles, Apj. 263, L1-L5, (1982). 

[12] G. R. Blumenthal, H. Pagels, J. R. Primack, Nature. 299, 37-38, (1982). 
[13] M. Tegmark. et. al., Astrophysics. J., 606, 702, |arXiv:astro-ph/0310725| 
[14] D. N. Spergel. et. al., Astrophysics. J. S., 170, 377, |arXiv:as tro-ph/0603449 
[15] J. R. Bond, A. S. Szalay, The Astrophysical Journal, 274:443-468, (1983). 



17 



[16] K. Sigurdson et al., Dark matter electric and magnetic dipole moments, Phys. Rev. D. 70, 083501 (2004). 
[17] S. D. McDermott et al., Turning off the Lights:How dark is darkt matter?. Phys. Rev. D. 83, 063509 (2011). 
[18] R. A. C. Croft, D. H. Weinberg, M. Pettini, L. Hernquist, and N. Katz, Astrophys. J. 520, 1 (1999). 
[19] Efstathiou et al., Nature. 348, 705-707, (1990). 
[20] Ostriker, Steinhardt, Nature. 377, 600-602, (2002). 

[21] P. J. E. Peebles, The Large scale structure of the Universe, Princeton University, Princeton, NJ, 1980). 

[22] H. Mo, F. C. Vanden Bosch, and S. White, Galaxy formation and evolution, Cambridge university, Cambridge, UK, (2010). 

[23] J. Jeans, Phil. Trans., 129, 44 (1902). 

[24] Astronomy and Cosmogony, Cambridge University Press, (1928). 

[25] E. Lifshitz, About gravitational stability of expanding world, Journal of Physics USSR. 10, 166 (1946). 

[26] E. M. Lifshitz, I. M. Khalatnikov, Advan. Phys. 12, 46, 185 (1963). 

[27] W. B. Bonnor, Monthly Notices, Roy. Astron. Soc, 117, 104 (1957). 

[28] M. Guyot, Ya. B. Zel'dovich, Astron. Astrophys. 9, 227-231 (1970). 

[29] P. Me'sza'ros, Astron. and Astrophys. 37, 225-228 (1974). 

[30] D. N. Spergel, J. S. Steinhardt, Phys. Rev. lett. 84, 3760, (2000). 

[31] V. Mukhanov, Gravitational Instability in Newtonian theory, Physical Foundations of Cosmology, (Cambridge University 
Press, 2005). 

[32] W. B. Bonnor, Monthly Notices, Roy. Astron. Soc, 117, 104 (1957). 

[33] S. Dodelson, The Boltzmann Equations, Modern Cosmology (Academic Press. Imprint of Elsevier, 2003). 
[34] R. K. Sachs, A. M. Wolfe, ApJ, 147-73S 

[35] W. HU, N. Sugiyama, The Astrophysical journal, 471, 542-570, (1996). 
[36] R. R. Caldwell, Elsevier. Phys. Lettr. B. 545, 23-29, (2002). 

[37] J. A. Frieman, M. S. Turner, D. Huterer, Annu. Rev. Astro. Astrophys. 46, 385-433, (2008). 



